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Abstract 

We numerically study the dynamical excitations in Bose-Einstein condensate (BEC) placed in 



periodic and quasi-periodic 2D optical lattice (OL). In case of the repulsive mean- field interaction 

Oh 

the BEC quantum tunnelling leads to a progressive soliton's splitting and generating of secondary 

•l-H 

solitons, which migrate to closest trapping potential minima. A nontrivial soliton dynamics appears 
when a series of 7r-pulses (phase kicks) are applied to the optical lattice. Such sudden perturbation 
produces a dynamic redistribution of the secondary solitons, leading to a formation of an artificial 
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solitonic superlattice. Different geometries of OL are analyzed. 
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I. INTRODUCTION 



Different types of solitons in Bose-Einstein condensates (BECs) have been intensively 
studied both theoretically and experimentally [If- [22]. Recently, solitons supported by weak 
attractive interactions between atoms were created in the condensate of 7 Li trapped in 
strongly elongated traps[23], [21], [2H]. Ref. [2E] detailes a set of experiments showing 
evidence of the existence of BEC of polaritons. Above a critical density, the authors ob- 
served massive occupation of the ground state developing from a polariton gas at thermal 
equilibrium. 

Also, a possibility to design ID and 2D quantum systems by "freezing out" one or two 
degrees of freedom by adding a ID or 2D optical lattice to the magnetic trap has been 
demonstrated experimentally [21], [28J . However the soliton's dynamics in two dimensions 
is much more involved and cannot be reduced to the one-dimensional case. 

It is well known that solitonic solutions appearing in integrable models, such as, for 
instance, ID spatial solitons in the cubic nonlinear Schrodinger equation (NSE)|2S], reveal 
particle-like behavior. In particular, they remain unperturbed after collisions, completely 
preserving their identities and form. However, there are no loosely bound stable solitons in 
2D (see e.g. [30]). 

The properties of solitons are well studied for a long optical lattice case (gap solitons). In 
the repulsive condensate gap solitons (GSs) can be created as a result of the interplay of the 
self-defocusing nonlinearity and periodic potential induced by an optical lattice (OL) - an 
interference pattern created by counterpropagating laser beams illuminating the condensate, 
[3T]-[35]. 

The dynamics of BEC in 2D optical lattice are currently a subject of intensive theoret- 
ical and experimental studies [36] -[H]. An interesting direction consists in using trapping 
quasiperiodic potentials for generation of several nonlinear excitations. The authors [12] 
have constructed some families of solitons for one- and two-dimensional Gross-Pitaevskii 
equations with a repulsive nonlinearity and a potential of quasicrystallic type (in the 2D 
case, the potential corresponds to a fivefold optical lattice). Some stable 2D and 3D solitons 
in the self-attractive Gross-Pitaevskii equation were found in [36] by using the variational 
approximation and direct simulations in real and imaginary time. Also, investigation of 
stability of localized vortices in an effectively 2D ("pancake-shaped") trapped BEC with 
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negative scattering length was performed in [39] . 

Usually, in experiments on solitons generation in BEC it is possible to obtain a unique 
soliton which is stabilized in gap zones of a periodic OL. However from the point of view 
of quantum-information technologies the possibility of generation of nonlinear excitations 
with complex structure, similarly to multisoliton state is more important. Several nonlinear 
regimes of interaction of such excitations have been studied in the literature. Nevertheless, 
the problem of experimental generation of multisolitonic states in BEC is still far from being 
completely solved. 

In this paper we explore an intermediate case when the size of BEC is about ten (or less) 
spatial periods of the optical lattice. In this situation the approximation of a perfect potential 
periodicity, V(x + L) = V(x) (necessary for the soliton stabilization) strictly speaking, is 
not valid, and a more realistic model should be constructed. In such configurations the 
dynamics of a single soliton are defined by the structure of the trapping potential in the 
vicinity of the soliton, rather by the far OL order. One can say that a soliton is subjected 
to a local spatially non-uniform nonlinear force field, which is defined by OL configuration, 
and by the excitation's amplitude. Such a field drastically changes the soliton dynamics: 
the soliton can tunnel into the closest minima (in the case of repulsive interaction), or turn 
to a self-compression dynamics (in the case of attractive interaction). 

We focus on the repulsive interaction case. Then, if the initial soliton position does not 
coincide with the minima of OL, the shape of such excitation will be split into fragments 
tending to flow into the closest minima of trapping potential, and the soliton dynamics 
essentially depends on a local spatial structure of OL. 

This observation suggests the exploration the soliton evolution in BEC when the spatial 
OL configuration is rapidly changed (which can be easily achieved experimentally). We 
have found that remarkable dynamics occur in situation when a soliton, steady in an initial 
OL configuration, becomes nonequilibrium in a new OL. This establishes the possibility to 
drastically increase the speed of forming the secondary solitons with such OL switching. 
Below we demonstrate that as a result of such a dynamics an artificial solitonic superlattice 
in BEC is formed. 

In this paper we analyze dynamics of a soliton placed in periodic or quasi-periodic optical 
lattices in 2D Bose-Einstein condensate (BEC) in the intermediate case (when the spatial 
scale of the condensate area is about 10 or less OL wavelengths), when the OL is suddenly 
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switched after applying periodical it phase shifts (phase kick). It will be shown that if 
the period of switching is sufficiently large (with respect to the tunneling time) secondary 
solitons are formed in OL minima. At later (successive) 7r -switching such secondary solitons 
pass through series of metastable states and migrate from the center to the periphery. In a 
suitable moment this process can be interrupted, which allows a free interaction of generated 
nonlinear excitations in agreement with the fundamental 2D nonlinear Schrodinger equation. 

The paper is organized as follows. In Section II we discuss basic equations for BEC 
atoms placed into OL. In Section III we present numerical studies of the soliton dynamics 
in the repulsive mean-field interaction for 2D case with periodical 7r-phase OL shifts. Var- 
ious geometries of OL, including periodic and quasiperiodic, are analyzed. We also discuss 
some peculiarities of the evolution in the attractive interaction case. In the last Section we 
summarize our results. 

II. BASIC EQUATIONS 

In order to trap a Bose-Einstein condensate in a (quasi)-periodic potential, it is sufficient 
to exploit the interference pattern created by two or more overlapping laser beams and the 
light force exerted on the condensate atoms. We mainly focus on the physical situation 
when the number of atoms is sufficiently large, so that atomic number fluctuations are 
negligible and the mean-field approximation can be applied. In this approach an anisotropic 
BEC cloud, loaded into a two-dimensional optical lattice potential V(r), is described by the 
Gross-Pitaevskii (GP) equation (T = 0) 

^ = [-vi + nr)+G|$| 2 ]$, (1) 

where r = {x,y}, Vj_ = {d x , d y }, $ is the condensate wave function, G = +1, —1 for 
repulsive and attractive interaction correspondingly. This equation is obtained by assuming 
a tight confinement in the direction perpendicular to the lattice ("pancake" trapping geom- 
etry) and the standard dimensionality reduction procedure (see e.g. [6], [L7]). According to 
Eq.Q, a particle with the (condensate) wave function $(r, t) evolves in the external poten- 
tial V(r) induced by OL plus the mean-field potential created by the remaining particles. 

It is made dimensionless by using the characteristic length h = pj ^8na s , energy E = 
g/p 2 = h/to, and time to = Hp 2 / g = mp 2 /4nhas, here and below m and Nq are the mass 
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and number of the trapped atoms respectively, g = 4irh 2 a s /m is the interaction strength, 
and a s is the s -wave scattering length. The wave function is made dimensionless as $ — ■> 
p$, where the factor p is specified by the normalization conditions N = j \p$>\ 2 dVo (i.e., 
no(r, r) = |p$(r,r)| 2 ), where Vo is the volume of the condensate. 

We consider a quasiperiodic time- dependent trapping potential V of the following form, 



N N 

V = V(x, y,r) = V c + eJ2 cos(k (n) r + tt0(t)) =V c + e(l - 26{r)) cos(k (n) r), (2) 

n=l n=l 

where kw = {cos(27r(n — 1) /N,sm(2n(n — 1)/N}, V c = const, and the time-dependence 
appears through the pulse function 6{r) = 1 for (2n + l)T t < r < (2n + 2)T t and 6(r) = 
otherwise, T t is a time period, n = 0,1, 2.... For stationary OL (T t — > oo) and iV = 5 such 
a potential case was studied in [12]. Some typical time independent potentials, V^x,?/, 0), 
are shown in Figs{T] (b),(c),(d) for the iV = 4, 5, 7 cases correspondingly. Such quasiperiodic 
optical lattices can be created in the 2D combination of N = 4, iV = 5 (Penrose 

tiling - a pattern of tiles, which completely cover an infinite plane in an aperiodic manner) 
or iV > 7 quasi-l-D sublattices with wave vectors k.^. The band-gap spectrum of 2D 
photonic crystals of the PT type has been studied in [33], [33], [35]. Recently the interesting 
properties of such lattices (e.g., Penrose lattices), have been discovered for quasiperiodic 
pinning arrays [36]). 

From now on we will concentrate on a repulsive interaction case, G = 1. First of all, we 
note that the generalized momentum density is not conserved because the trapping potential 
V(x, y) induces a spatial inhomogeneity in a system [49] . Let us estimate the force acting 
on the condensate particles. It is well known that the GP equation (fl]) can be generated 
from the Lagrangian density 

L = ^[$*$ T - - | V ± $| 2 - V(r) |$| 2 - - |$| 4 , (3) 
and the corresponding full condensate energy has the form 

E = J dr[|V ± $| 2 + V(r) |$| 2 + - \§\ 4 }=const. (4) 

Observe that the constant part of the trapping potential, V c , can be removed from the 
Lagrangian ^ by renormalizing the wave function as follows, 




FIG. 1: (a) Init Gaussian package. (b),(c),(d) two-dimensional distribution of the trapping poten- 
tial V(x, y, 0) induced by OL, Eq.Q for cases (b) N = 4; (c) N = 5; (d) N = 7. 



$ — > $exp(-iV c r). (5) 

Thus, without loss of generality we suppose that the potential V(x,y) in ([I]) and ([3]) can 
be renormalized such that V min (x,y) = and V(x,y) > 0. Note that in the case of such 
renormalization the energy Q is shifted by V C N , which allows us to associate V c with the 
chemical potential \x = dE/dN . 

The Lagrangian (|3| generates the following force field acting on the soliton, 

F ^ = 6¥ = d¥- v d^ = Fv + Fs + F ^ (6) 

where 

F v = -V(x, y, r)$, F N = -G |$| 2 $, F 5 = vi$. (7) 



6 



One can observe that Fy = in the vicinity of V m i n = and therefore in the potential minima 
the soliton evolution is governed by the standard nonlinear Schrodinger model dynamics, i.e. 
only subjected to the forces F s and F N . This, in particular, leads to a very low tunneling 
rate even for an intermediate depth of OL. Thus, as we have numerically corroborated, a 
sufficiently narrow soliton placed in one of the potential minima is practically a steady-state. 

It is clear that a deviation from V m i n = generates a nonzero Fy, which essentially 
modifies the soliton behavior. The effect of Fy becomes predominant in the vicinity of the 
maximum of the potential. Thus, after applying a phase kick (when the minima of the 
potential are suddenly converted into maxima), the initial "steady-state" soliton starts to 
rapidly evolve in the direction of the closest minima of the new potential. Besides, in the 
vicinity of the soliton maximum (where $ is real and positive), the forces Fn and Fg are all 
negative. These forces stretch and deform the initial packet leading to the soliton splitting. 
It is worth noting that due to the energy conservation Q, the decay of some part of the 
BEC soliton amplitude |$| must be compensated by sharp gradients of some other parts. 



III. NUMERICAL RESULTS 
A. Repulsive interaction 

Previous arguments have only heuristic meaning and do not provide a quantitative de- 
scription of the soliton dynamics in the 2D case, especially for time-dependent OL potential. 
The deeper insight into the BEC dynamics in a general OL depth case can be reached only 
by numerical methods. We have investigated two-dimensional dynamics for various ini- 
tial conditions and different parameters of the trapping potential. The results are shown 
in Figsj2]j6j For our numerical experiments we used the combination of known splitting 
methods[50j generalized to the nonlinear case. We have numerically solved Eq.Q in do- 
main 200 x 200 with zero boundary conditions for times r < 200. Greater temporal intervals 
normally require a spatial grid with greater size. Our calculations were done in the complex 
$ plane to observe the evolution of both amplitude and phase of the soliton. Both the norm 
of the wave function and energy Q were conserved with good accuracy (less then 1%). We 
have studied the evolution of the initial Gaussian packet, FigjlVa), 
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$(x, y, 0) = $ exp[-q x (x - x ) 2 -q v (y- y f) (8) 

for different factors q XjV , and the amplitude $o> which is a function of the total number of 
trapped atoms N . We have numerically analyzed the packet ^ evolution initially centered 
at different points of the plane (x ,y ), corresponding either to minima (y min ) and maxima 
(y max ) of the potential. 

Let us start with a periodic trapping potential Eq.Q, corresponding to iV = 4 (see Figjl] 
(b)) and for G — 1. We have used the parameters: q x = q y = 0.004, e = 0.7, $o = 1, and 
the initial packet was placed at the central minimum, FigJTJa). In the case of a stationary 
potential the initial packet evolves in a well known way by splitting into fragments and 
slowly tunneling into the closest minima of the optical lattice. Nontrivial dynamics arise 
after sudden (short with respect to the scale to) shift by 7r of the optical lattice phase, when 
OL minima are switched with their maxima. Being placed in these new maxima the soliton 
flow out into new closest minima and after a certain transition time the secondary solitons 
are generated in those minima. In our simulations the period of the phase shift, T t = 50, is 
chosen larger than the above mentioned transition time (proportional to the characteristic 
time of quantum tunneling), so that the details of the initial state are washed out and the 
secondary "solitons" acquire well pronounced shape. 

These dynamics are shown in Fig|2j In Figj2] (a) the soliton state at time r = 40 before the 
first phase switching is shown. One can observe the beginning of slow condensate tunneling 
into closest potential minima. In Fig|2] (b) the condensate state at r = 80 (after first phase 
switching) is presented. One can appreciate an essential reduction and widening of pulse 
amplitude. We observe that such dynamics become significantly faster with respect to the 
stationary case. At r = 120 (after second phase switching), Figj2] (c), the state acquires a 
more complex shape with quite pronounced secondary solitonic peaks at the corners. The 
final state at r = 160 (after third phase switching) is shown in Figj2] (d). We observe already 
a well formed artificial solitonic superlattice with the peaks located in the trapping potential 
minima. It is clear that the details of this new soliton structure depend mainly on the local 
vicinity of OL. 

Due to a recent interest to quasiperiodic structures, we also analyze the BEC evolution 
in an optical structure with a lower symmetry trapping potential. Such dynamic is shown 
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FIG. 2: Evolution of BEC soliton for period 7r-phase shift Tt = 50 and different times: (a) r = 40; 
(b) r = 80; (c) r = 120; (d) r = 160. See details in text. 

in Fig|3]for potential (|2]) with N = 5 and in Fig|4]for potential ^ with N = 7 (see Figjl] (c) 
and (d) accordingly). The case N = 5 is essentially a two-dimensional Penrose-tiled lattice, 
where the tiles are two kinds of rhombus: a thin tile (with vertex angles of 36° and 144°) and a 
fat tile (72° and 108°) [H]. We observe from Figj3]and Figjijthat for quasiperiodic potentials 
the system evolves in a less symmetric way. The condensate is distributed on closed smooth 
cells without strongly pronounced solitonic peaks. Such cells appear more smooth in the 
5-fold case, see Fig|3j This effect is similar to 2D Penrose lattices in superconductors, where 
the pinning of vortices is related to matching conditions between the vortex lattice and the 
quasiperiodic lattice of pinning centers 06]. We notice again that the final structure of a 
condensate still has the form of localized states distributed on the minima of the potential. 

It is representative to illustrate the state evolution as the time behavior of the amplitude 
|$| of the maximal peak, as shown in Figj5^a). It can be observed that the amplitude after 
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(a) ft) 




FIG. 3: The same as in Figj2] except N = 5 in trapping potential V(x, y, r) Eq.Q. 

some period of growing (which is accompanied by the stretching of the initial distribution 
due to non-linear effects) starts to diminish and reaches its minimum (quasistationary value) 
at r = 130, although the splitting is still continuing. 

To evaluate quantitatively the behavior of BEC shape we have studied the behavior of 
the soliton "center of mass" R(r) defined as 




where ro indicates the initial soliton position. The dynamics of R(t) = |R(t)| are shown 
in Figj^b). One can observe that R(r) grows in a very different way for periodic and 
quasiperiodic potentials, which is a reflection of the different local structure of OL. On the 
other hand, while the minima are progressively filling up with the condensate the quantity 
R(t) tends to saturation in both cases. 
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FIG. 4: The same as in Figj2] except N = 7 in trapping potential V(x, y, r) Eq.Q. 
B. Attractive interaction. 

Although the main goal of this paper is studying the repulsive interaction, for com- 
pleteness, we also discuss the dynamics of the initial Gaussian pulse ^ in the attractive 
interaction case (G = —1). Let us recall that for the 2D nonlinear Schrodinger model (with- 
out potential) the dynamics of the average radius of an excitation are (see e.g. Chap. VIII in 
Ref-EHU) (R (r)) 2 = J p 2 |$| 2 d 2 r ~ (E/N )t 2 + R 2 (see also^T] and references therein). So, 
if the energy Q is negative, E < 0, then such an excitation undergoes a self-compressing. 
For the initial Gaussian ^ this yields the inequality $g > g Xiy , which is satisfied in our case. 
However a similar analytical condition for 2D BEC with OL is not established yet in the 
literature. Therefore we have studied the soliton evolution for the 2D attractive interaction 
numerically for the initial state FigjT] (a) and different OL geometries (Figjl](b,c,d)), similar 
to the repulsive case. 
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FIG. 5: Time dynamics in case of repulsive interaction (G = 1). (a) |$ max (T)| / |^max(0)|, and (b) 
mass center R(t), see ([9]). Solid line, dash line and dash-dots line correspond to N = 4, N = 5 
and N = 7 cases accordingly for trapping potential V(x, y, r) in Eq.Q. Here r = n (n is the time 
step index). 



As it can be appreciated from Fig|6] the soliton dynamics completely change in the at- 
tractive case G — — 1. It is well known that for the 2D attractive interaction the collapse 
dynamics takes place and there are no loosely bound stable states. We have found that 
for our initial state the collapse dynamics take place and there are no loosely bound stable 
states in two dimensions. Nevertheless, it turns out that such instabilities can be exploited 
in order to prepare BEC solitons[T7]. We have investigated such a dynamics for various 
geometries of the trapping potential. In contrast to the repulsive interaction case we found 
that the self-compressing behavior is practically insensible to changing OL symmetry. We 
observe from Fig|6]that just a single pulse is generated from the initial packet, which evolves 
without any splitting or visible deformation. Such self-compression leads to a singularity 
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FIG. 6: Self-compressing of initial Gaussian package at attractive interaction G = — 1 and the 
period of vr-phase shift T t = 8. (a) r = 15; (b) 20; (c) 25; (d) 30. 

and therefore, the local structure of OL in the vicinity of the soliton does not notably affect 
its evolution. Obviously, the tunneling processes quickly become practically negligible in 
this case. We have also found that small shifts of the soliton center from the minimum of 
the potential lead to the same dynamics, and in particular, do not change the single-peak 
structure of the final state. Generating of two and more peaks as well as stabilization of 
self-compressing can be reached via a specific initial state and/or OL geometry. However in 
general this question remains open and will be studied elsewhere. 



IV. CONCLUSION 

We have studied the dynamics of soliton in Bose-Einstein condensate (BEC) placed in 
periodic or quasi-periodic 2D optical lattices (OL). In the case of the repulsive mean-field 
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interaction the BEC tunneling leads to progressive deformation and splitting of solitons. It 
results in the generation of secondary solitons migrating to closest trapping potential minima. 
We found that nontrivial dynamics arise when a series of 7r-phase shifts (phase kicks) are 
applied to the OL, so that the minima and maxima of OL are periodically exchanged. This 
results in a progressive migration of the solitons from the initial position to the periphery 
with the forming of the artificial solitonic superlattice. We have found that the geometry of 
forming solitonic superlattice is rather sensitive to the geometry of the trapping potential. 

Such an effect suggests an effective and simple method for splitting of the initial BEC 
soliton and generation of multisolitonic nonlinear excitations in the experiment, and thus, 
can be, in principle, applied in quantum information technology, since it allows a creation 
and manipulation of solitons as bits (qubits) in BEC quantum structures. 
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